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INTRODUCTION 


For  the  purpose  of  energy  storage  we  seek  a  molecule  such  as  N2  or  CO  which  has 
one  or  more  we^y  bound,  yet  possibly  stabilizable,  upper  states  and  a  very  deeply  bound 
ground  state.  Figure  1  illustrates  the  concept.  It  shows,  for  example,  that  one  might  be 
able  to  extract  alraut  1 1  eV  for  each  transition  of  the  weaJdy  bound  excited  state  of  CO  to 
the  ground  state. 


Figure  1.  Potential  Energy  Release  frmn  a  High  .Spin  , 

High  spin  state  molecules  corresponding  to  the  interaction  of  ground  state  atoms  are 
particulariy  likely  candidates  for  the  upper  state  since  they  typically  will  exhibit  mainly 
electrostatic,  induction,  attd  dispersion  type  van  der  Waals  interactions  at  long-range  but  tK> 
short-range  "chemical"  binding.  Furtheraore,  such  high  spin  states  will  have  no  dipole- 
allowed  transitions  to  the  ground  state  and  thus  are  potentially  long-lived  states.  Recently, 

Hzpptx  and  coworkers  (Ref.  1)  have  argued  that  N2  (and  by  implication  CO) 

has  spin-spin  decay  nnodes  that  cause  it  to  have  an  unacceptably  short  lifetime.  Since  these 
decay  processes  apparently  depend  on  R~^  where  R  is  an  interfiragment  distance,  molecular 
size  is  an  important  parameter.  If  we  can  increase  the  effective  size  of  the  ^1.*  CO  by 
surrounding  it  with  a  cluster  of  H2  molecules  (bound  to  CO  by  van  der  Waals  forces), 
however,  we  may  be  able  to  increase  the  high-spin  lifetime  enough  so  that  the  cluster 
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becomes  an  interesting  high  energy  candidate.  Thus,  we  have  launched  a  detailed  study  of 
the  quintet  states  of  the  CO  molet^e  that  has  the  following  main  conqx}nents; 

(1)  A  preliminaiy  survey  of  the  long-range  interactions  of  ground  state  C  and  O  atoms 
by  means  of  a  semi-enqjtrical  m^l; 

(2)  Detailed  high-quality  ab  initio  computations  of  the  potential  energy  curves  for  the 
quintet  states  of  CO  which  correspond  asymptotically  to  ground  state  atoms; 

(3)  Detailed  spin-orbit  and  spin-spin  studies  to  determine  the  decay  modes  of  the  high- 
spin  CO,  and  to  learn  whether  the  high-spin  lifetime  can  be  increased  by 
surrounding  CO  with  a  cluster  of  H2  or  He. 

For  the  sake  of  clarity  and  conqrleteness,  this  report  will  include  portions  of  the 
research  which  had  been  performed  un^r  contract  F0461  l-86-K-0<)75/P00001  let  to  the 
State  University  of  New  York.  Some  of  the  results  of  our  researdi  have  been  prepared  in 
form  appropriate  for  publication  in  Chemical  Physics.  A  semi-final  idraft  of  the  journal 
article  appears  as  an  appendix  to  this  report.  In  the  first  section  of  this  report  under 
"Technictd  Discussion"  I  present  a  summary  of  die  journal  article  which  also  serves  as  a 
prologue  to  the  nh  initio  investigations. 

The  next  sectirm  under  "Technical  Discussion”  comprises  a  report  of  our  ab  initio 
quantum  mechanical  investigations  of  the  quintet  states  of  CO  which  correspmid  to  the 
interaction  of  ground  state  carbcm  and  oxygoi  atoms.  This  is  cmnponent  (2)  mentioned 
above.  First  we  introduce  die  problem  of  the  high-^>in  states  of  CO  and  give  a  review  of 
the  very  scanQr  badcground  literature.  Then  we  d^elcp  the  basis  set  which  is  optimized  to 
reproduce  relevant  energies  and  electric  mult^le  polarizabilities  of  die  carbon  and  oxygen 
atoms.  Next  we  descrite  the  ab  initio  conqnitations  of  the  potential  energy  corves  for  the 

1,2  1  ^A,  1  ^Z~,  and  1,2  ^11  states  of  CO,  and  finally  we  discuss  die  results.  The 

summary  indicates  the  direction  of  additiorud  study  which  is  required  to  estimate  the 
lifetime(s)  of  these  quintet  states.  That  will  be  con^ndnt  (3),  mentioned  above,  of  the 
study  of  CO. 
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TECHNICAL  DISCUSSION 


The  Long-Range  Interaction  of  Ground  State  Carbon  and  Oxygen  Atoms  in 
the  Zero-Overlap  Approximation 

For  our  present  purposes  a  reasonable  review  of  the  literature  of  long-range 
interatomic  interactions  is  given  in  a  few  key  references  (Ref.  2-S).  Suffice  it  to  say  that 
die  great  bulk  of  the  treatments  of  long-range  interactions  involve  either  (a)  an  empirical 
qiproach  where  the  form  of  the  interaction  potential  is  assumed  (perhrps  guided  by  t^ry) 
and  whose  parameters  are  fixed  by  fitting  to  experimental  data  of  various  kinds,  or  else  (b) 
a  strictly  theoretical  qiproach  where  the  energy  of  interaction  is  calculated  either 
variatiotiably  or  using  perturbation  theory.  Typic^y,  the  latter  qiproach  has  been  limited 
to  the  interactimi  of  ground  state  species.  Recendy,  Bussery  and  Aubert-Frecon  (Ref.  6-8) 
have  developed  a  semi-enqiitical  method  for  estimating  the  Interaction  energies  of  atoms  in 
the  "zeto-overlq)''  regime  of  interatomic  separations.  Hieir  scheme  (ixbich  we  will  call  the 
BAF  scheme)  overccnnes  certain  limitations  of  nuuiy  theoretical  approaches,  which  allows 
it  to  be  used  to  describe  excited  states.  The  BAF  scheme  follows  the  usud  perturbation 
tfaecvy  develqm^  (Ref.  2)  to  the  point  where  one  needs  to  evaluate  the  cooqwnents  of  the 
energy  expression  in  the  multipole  expansion.  They  recognize  that  they  can  use 
experimental  transition  moment  data  and  values  of  experimental  atomic  energy  levels 
instead  of  calculated  values.  When  necessary  they  calt^ate  excited  atomic  energies  by 
using  Whittaker  functions  to  represent  the  wavefiinction.  BAF  have  applied  their  scheme 
to  various  halogen  or  alkali  diatomic  molecules  with  quite  notable  success  as  judged  by 
agreement  with  high  quality  spectroscopic  observations  or  high  quality  ab  initio  results. 

Bussery  and  coworkers  (Ref.  9)  have  qiplied  the  BAF  scheme  to  the  interaction  of 
ground  state  C(^Pj)  and  0(^Pj)  atoms.  They  find  that  the  long-range  interaction  energy 
may  be  represoked  by  the  inult^le  eqiression 

E(R)  -  C3  R-5  +  Qi  R-<>  +  C8  R-*  +  Cio  R-^O  +  ...  (1) 

The  C5R~^  term  represents  die  energy  contribution  doe  to  the  electrostatic  interaction  of  the 
permanent  quadrupole  in  C  with  d^  in  O;  that  term  may  be  either  positive,  negative,  or 
zero  depending  on  the  coupling  of  the  atomic  angular  momenta  to  produce  a  particular 
molecular  state.  The  remaining  terms  in  Eq.  (1)  represent  the  various  components  of  the 
dispersion  interaction  energy  which  corresponds  to  the  interaction  of  mutually  induced 
mult^les.  For  examine,  die  term  represents  the  interaction  of  an  induced  d^le  on 
C  wifo  an  induced  dipole  on  O.  The  magnitude  of  Q  dqiends  on  the  particular  state  which 
results  from  the  interaction  of  the  atoms,  but  its  sign  is  intrinsically  negative  (thus,  the 
interaction  is  always  attractive)  ordy  so  long  as  both  interacting  atoms  are  in  their  ground 
electronic  states.  The  higher  order  terms  correspond  to  the  interactions  involving  higher 
order  induced  multipoles. 

Bussery  and  coworkers  (Ref.  9)  argued  that  their  values  of  C5  are  relatively  reliable 
but  because  the  experimental  values  of  transition  moments  and  energy  levels  needed  to 
calculate  C6  in  dw  BAF  scheme  ate  so  uncertain,  diey  can  only  obtain  erode  estimates  for 
C5.  They  did  not  even  calculate  values  of  Cg  or  Cio  according  to  the  BAF  scheme  since 
die  results  for  the  higher  coefficients  were  likely  to  be  even  more  uncertairu  They  did  find, 
however,  that  the  €5  value  was  roughly  the  same  for  all  molecular  states  which  correspond 
to  the  interaction  of  ground  state  C  and  O  atoms. 
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Since  it  is  clear  that  cmly  a  cnide  estimate  of  the  zero  overly  interaction  energy  may 
be  obtained  presently  in  any  case,  we  have  obtained  some  rough  estimates  of  the  C^,  Cg, 
and  Cio  coefficients  by  using  the  highly  approximate  expressions  13.3.36  -  13.3.38  of 
Ref.  1.  They  are 


Ce  =  (3/2)  Ic  lo  <xcao/(Ic  +  lo). 

(2) 

♦  Jio  +  icJ 

(3) 

Cio  =  (^^j(lcloocao)2  /ac  +  lo). 

(4) 

where  e  is  the  charge  of  an  electron,  and  where  we  have  substituted  in  the  corresponding 
expressi<»is  given  in  Ref.  1  the  values  (a)  the  first  ionization  potential  of  the  atoms  for  the 

term  hu  and  (b)  die  average  value  of  the  diptde  polarizability  calculated  in  the  section  of  this 

report  for  the  a  terms.  We  compare  in  Table  1  the  values  of  the  Cn  of  Eq.  (1)  given  by 
Bussery  and  cowoikers  with  drose  calculated  frmn  Eqs.  (2-4). 


TABLE  1.  Long-range  Coefficients  of  Eq.  (1)  in  a.u. 


Source 

C5 

C6 

C8 

Cio 

presem 

11.2 

19.7 

317 

2950 

Bussery  and  cowotkers 
(Ref.  9) 

11.2 

13.6 

352 

In  any  case  die  long-range 

analysis  suggests 

that  all  the  quintet 

states  of  CO  with 

the  possible  exception  of  the  2  ^11  state  will  be  attractive  in  the  long  range.  Each  attraction 
will  be  balanced  by  the  exchange  repulsion  at  some  intermediate  separation  which 
consequently  gives  a  minimum  in  each  curve.  At  sufficiently  small  separations  the  quintet 
potential  curves  will  be  repulsive.  No  matter  which  estimate  of  the  terms  one  chooses 

to  use,  the  long-range  analysis  suggests  that  the  1  will  be  by  far  the  most  deeply  bound 
of  the  quintets  corresponding  to  C(^Pj)  CK^Pj).  The  interaction  curves  will  lie  in  the 
order  1  «  1  <  2  1  1  5n,  <  2  at  sufficiently  large  separations.  In  the 

next  section  we  examine  these  predictions  by  means  of  more  sophisticated  theory  in  order 
to  learn  whether  any  of  the  quintets  is  sufficiently  deeply  bound  to  be  an  interesting  energy 
storage  candidate. 


Ab  InUio  Computation  of  the  Interaction  Energy  Curves  for  the  1,2 
1  1  and  1,2  *11  States  of  CO 


I  have  established  the  motivation  for  investigating  the  quintet  states  of  CO,  and  in 
the  last  section  I  gave  a  crude  analysis  to  show  that,  with  the  likely  «iception  of  the  2  *n 
state,  the  quintets  of  CO  are  likely  to  exhibit  long-range  (or  van  der  Waids)  binding  and. 
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thus,  are  possible  energy  storage  candidates.  Unfortunately,  there  is  presently  no 
experimental  spectroscopic  data  with  which  to  compare  these  crude  predictions.  Such  a 
situation  is  not  unique  to  CO.  Although  fluorescence  from  excited  sodiiun  vapor  was 
described  by  Wood  (Ref.  10)  about  80  years  ago,  and  a  few  of  the  lower  lying  transitions 
in  the  molecular  singlet  manifold  were  mapped  out  spectroscopically  about  50  years  ago  by 
Locnnis  and  Kusch  ^ef.  1 1),  it  was  quantum  dieoretical  computations  in  the  last  decrde  or 
so  (Ref.  12)  u^ch  first  gave  reliable  informatirni  on  molecule  states  not  readily  accessible 

fixmi  the  ground  X  state  of  Na2.  (And  Na2  is  a  much  simpler  molecule  than  CO!)  I 

am  aware  of  only  the  pioneering  quantum  mechanical  computations  of  O'Neil  and  Schaefer 
[ONS]  (Ref.  13)  that  even  discuss  the  quintet  states  of  CO.  Unfortunately,  their  work 
extends  only  to  an  intemuclear  separation  of  Sao,  and  die  Bussery  and  coworkers  estimates 
(Ref.  9)  are  valid  at  separations  no  smaller  than  7  Sq.  The  ONS  results,  which  ate  obtained 
from  a  full  configuration  interaction  treatment  based  on  a  minimal  basis  set,  suggest  that 

(Hie  or  both  of  die  1  and  1  ^11  states  are  s<Hnewhat  more  deeply  bound  than  one  might 
expect  on  the  basis  of  the  crude  lorig-range  analysis  discussed  in  the  last  section. 
Consetpiendy,  it  qipears  that  ihore  detailed  computations  and  analysis  are  in  order. 

Basis  Set  Development 

Recall  that  it  is  our  goal  to  obtain  accurate  potential  energy  curves  for  the  (piintet 
states  of  CO  which  correspond  asymptotically  to  die  interactkHi  of  C(3p)  and  0(3p).  Since 
we  presume  these  molecular  states  to  be  weakly  bound  van  der  Wa^  states,  it  is  doubly 
inqiottant  that  the  basis  sets  for  these  calculations  be  carefully  chosen. 

For  both  atoms  we  began  with  the  15  function  Slater  basis  sets  of  McLean  and  Liu 
[ML]  (Ref.  14)  which  are  given  in  Table  2.  The  starting  bases  were  augmented  by  a 
Rydberg  3s  and  a  RydRierg  3p  functitm  (so  that  all  the  molecular  calculations  will  have 
orbitals  with  scxne  Rydberg  character).  We  discuss  the  optimization  procedure  later.  We 
further  augmented  ari^or  (dunged  the  basis  to  include  polarizing  fuiK^ons.  The  latter  ate 
necessary  to  account  for  the  polarization  of  an  atom's  charge  distribution  when  another 
at(Hn  or  a  molecule  is  brought  near  it.  These  polarizati(Hi  effects  are  the  driving  forces 
behind  the  formation  of  van  der  Waals  interactions.  In  all  cases,  we  re-optimized  ^  most 
difiiise  3d  functions  of  the  origirul  ML  bases  for  the  dq>ole  polarizability,  we  optimized  an 
added  4f  function  for  die  quatlrapole  polarizability,  and  we  optimized  an  added  5g  function 
for  the  octupole  polarizability.  I^uthermore,  we  performed  a  test  for  the  carbon  atom  to  see 
if  adding  a  p  furiction  inqnoved  the  (]uadrupole  polarizability.  It  did  not.  In  a  similar  test 
we  added  a  d  function  to  see  if  the  octupole  polarizability  was  changed  significantly  by  its 
additkm.  The  octupole  polarizability  of  C  was  changed  by  about  3%  (frcnn  699  ao^  to  719 
tifP)  by  adding  a  3(1  fun^on.  This  was  not  deemed  a  large  enough  change  to  include  it  in 
an  already  very  large  basis.  However,  it  indicates  that  the  error  in  the  octupole 
polarizability  is  relatively  small.  The  additions  and  changes  to  the  ML  basis  are  given  in 
Table  2. 

A  few  words  about  the  optimization  procedure  are  in  order.  We  used  the 
ALCHEMY  (Ref.  15)  system's  Slater  integrals  and  SCF  programs  to  optimize  (i.e., 
minimize)  die  energy  of  Rydberg  states  of  the  atoms.  For  cafoon,  this  was  done  by 
minimizing  the  energy  with  respect  to  the  added  functions  of  the  triplet  P  states  with  the 

following  electron  configurations:  Iso^  2pa  3so^  2pn  (to  obtain  the  3s  functions)  and 

Iso^  Iso^  2so^  2pa  3pic  (to  obtain  the  3p  functions).  It  should  be  noted  that  these  are  not 
the  lowest  energy  states  of  carbtni  containing  3s  and  3p  orbitals.  However,  in  a  single 
configurati(m  calculation  the  2s(2p)  orbitds  have  to  be  removed  to  prevent  orbital 
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degeneracies  fiom  occurring.  For  oxygen,  we  added  the  3s  and  3p  functions  and 
optimized  Hk  SCF  energy  of  ^  quintet  S  and  P  states. 


TABLE  2(a).  Modifications  of  the  McLean  and  Liu  (Ref.  14)  6s,  4p,  3d,  2f  Slater  Basis 
Sets  for  C  and  O  to  Add  Rydberg  Character  ai^  IiTq)rove  Multipole 
Polarizabilities. 


scq  # 

n 

1 

C 

carbon 

values 

oxygen 

1 

1 

0 

9.4826 

13.7574 

2 

1 

0 

5.4360 

7.6141 

3 

2 

0 

4.2010 

5.8660 

4 

2 

0 

2.6844 

4.3120 

5 

2 

,  0 

1.5243 

2.4802 

6 

2 

0 

1.0575 

1.6982 

8 

2 

1 

6.5100 

7.5648 

9 

2 

1 

2.6005 

3.4499 

10 

2 

1 

1.4436 

1.8173 

11 

2 

1 

0.9807 

1.1439 

13 

3 

2 

3.6407 

4.8299 

14 

3 

2 

2.0211 

2.5442 

15 

3 

2 

1.3730 

1.6015 

16 

4 

3 

2.5985 

3.2711 

17 

4 

3 

1.7653 

2.0590 

TABLE  2(b). 

Rydberg  and  Polarizing  Functions.  The  sequence  numbers  indicate 
where  diese  fimctions  belong  in  the  basis  set  of  Table  2(a). 

7 

3 

0 

0.60186 

0.70536 

12 

3 

1 

0.48707 

0.54074 

15* 

3 

2 

0.87500 

1.24837 

18 

4 

3 

0.92500 

1.18750 

19 

5 

4 

0.94375 

1.16250 

*  This  replaces  basis  tanctkm  #15  from  the  original  McLean  and  Lin  basis. 


To  obtain  the  dipole,  quadiupole,  and  octupole  polarizabilities,  we  used  the  ATMBIS 
program  of  A.C.  Wahl  as  amend^  by  W.  J.  Stevens  (Ref.  16)  to  include  the  finite  field 
approximation  (Ref.  17).  The  added  (or  changed)  functions  were  diosen  to  maximize  the 
energy  difference  between  the  field  "on"  arxl  the  field  "ofT  cases.  The  maximization 
process  always  took  place  with  the  function  to  be  added  present  in  either  the  sigma  space 
(denoted  by  ^  =  0)  or  else  only  in  the  pi  orbital  space  (denoted  by  =  1)  (except,  of 
course,  for  the  optimization  of  the  3s).  We  determined  the  final  values  of  the 
polarizabilities  with  the  optimized  functions  present  in  all  the  qrpropriate  orbital  spaces. 
The  polarizabilities  and  atomic  energies  we  obtained  are  given  in  Table  3.  There  we  show 
that  both  the  sigma  (M^  »<))  and  pi  (^  »!)  projections  of  the  polarizidnlities  and 
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T  ABLE  3.  D^lc,  Quadiupole  and  Ocmpole  Polarizabilities  of  C  and  O  in  Atomic  Units. 


ML  Basis 

hSSSDl 

Mover  and  coworicers  a.b 

Od 

II 

o 

— 

10.102 

lO.lOa 

Od 

1 

12.362 

12.990 

13.05a 

Od 

— 

12.021 

12.07a 

Yd 

— 

2.888 

2.95a 

Oq 

o 

II 

41.60 

41.73b 

Oq 

1 

46,49 

61.80 

61.28b 

— 

55.07 

54.76b 

Yq 

— 

20.20 

19.55b 

Oo 

=0 

— 

449.57 

Oo 

1 

253.18 

697.73 

— 

Oo 

— 

615.01 

— 

Yo 

— 

248.16 

— 

Od 

Mf  *0 

— 

5.005 

5.14a 

Od 

1 

4.23 

4.540 

4.58a 

— 

4.695 

4.77a 

Yd 

- 

— 

-0.465 

-0.56a 

«q 

o 

II 

S' 

20.32 

20.6lb 

«q 

1 

11.41 

14.74 

15.05b 

Oq 

— 

16.60 

16.90b 

Yq 

— 

-5.58 

-5.56b 

Oo 

*0 

— 

141.542 

— 

Oo 

1 

37.33 

98.352 

— 

Oo 

— 

112.749 

— 

Yo 

— 

-43.190 

— 

a  Ref.  18 
b  Ref.  19 
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consequently  the  average  value  (a)  and  the  anistropy  (y)  agree  nicely  with  the  carefully 
optimized  values  of  Meyer  and  coworkers  (Ref.  18, 19).  The  energies  we  calculated  with 
these  basis  sets  are  given  in  Table  4.  There  we  show  that  the  atomic  excitation  energy 

values  (AE)  which  we  calculate  agree  reasonably  well  with  experiment  (Ref.  20). 

TABLE  4.  Atomic  Energies  (in  Hartree  Atomic  Units)  and  Excitation  EiMrgies  (in  eV) 
for  C  and  O  Atoms. 


AE 


Atom 

State 

-E(SCF) 

-ELCD 

Cl  Calc 

Expt. 
(Ref.  20) 

C 

3p 

37.6886303 

37.7800911 

Id 

0.6315801 

0.7316137 

1.32 

1.26 

is 

0.5770546 

0.6819760 

2.67 

2.68 

0 

3P 

74.8102776 

74.9767572 

ID 

0.7299571 

0.9032159 

2.00 

1.97 

is 

0.6507112 

0.8221247 

4.21 

4.19 

ralmlarinn 

Eidier  a  multi-configuration  self  consistent  field  (MCSCF)  or  a  configuration 
interaction  (Q)  calculation  is  needed  to  obtain  relatively  reliable  potential  energy  corves  for 
the  weakly  interacting  systems  of  interest  here.  At  large  intemuclear  separations  the 
molecular  oibitals  are  essentially  atomic  in  nature.  ITius,  we  can  represent  the  electronic 

2  2  2  2  '4  2 

structure  by  the  configuration  Isq  Is^  2sq  2s^  2pq  2p^,  where  we  have  written  the 

atomic  orbitals  in  the  ascending  order  of  the  approximate  energies  of  the  atomic  orbitals. 
Cileariy,  we  must  correlate  all  six  of  the  valence  level  p  electrons,  and  perhaps  also  correlate 
the  cartwn  2s  electrons  (for  a  total  of  eight).  We  deemed  it  to  be  too  expensive  to  correlate 
the  oxygen  2s  electrons  as  well  for  a  to^  of  laa  correlated  electrons. 

Since  the  treatment  of  the  1  ^A  and  1  states  is  quite  straightforward,  let  us 
discuss  those  calculations  first.  We  use  die  ouqnit  vectors  firom  an  SC!F  calculation  for  the 
1  ^1'*'  state  as  input  to  the  valence  MCSCT  calculation  with  six  active  electrons  distributed 
among  two  sigma  and  two  pi  orbitals.  We  call  this  a  {6c~f2a  -f  2k)  MC  computation.  The 

MC  computation  for  the  state  comprised  four  configuration  state  functions  (CSF).  We 
used  the  MC  output  as  the  input  vectors  for  a  second  order  CT  (SCXH)  calculation  which 
comorised  all  single  and  double  excitations  of  eight  electrons  fixnn  three  a  and  two  n 
valence  orbitals;  we  call  these  a  (8e~/3o  -1-  2k)  CT  computation.  It  comprised  386389 
CSF.  We  also  checked  to  see  how  the  8e~  CT  results  were  affected  by  basing  them  on 
vectors  obtained  from  an  (8e~/3o  +  2n)  MC  calculation  (which  conqxrised  21  CSF). 


The  (8er/3<y  +  2n)  Cl  binding  energies  for  the  state  were  insensitive  to  whether 
we  used  the  (6e~/2o  +  27c)  or  the  (8e“/3a  +  2n)  MC  vectors  as  input.  The  results  listed  in 
Table  5  result  from  the  (8e“/3o  +  27c)  SOCI  based  on  (8e~/3cy  +  27i)  MC  vectors.  The 
computations  for  the  state  were  quite  similar  to  those  for  the  state.  The  %  results 
given  in  Table  5  correspond  to  an  (8e“/3a  +  2ji)  SOCI  computation  comprising  540,163 

CSF  and  which  itself  is  based  on  an  (8e~/3a  +  2it)  MC  computation  comprising  21  CSF. 
(The  SOCI  computation  step  for  each  point  required  about  6  h  of  CPU  time  on  an  IBM 
3090  computer.). 


T  ABLE  5.  The  Potential  Energies,  in  Units  of  -Efa,  of  Low  Lying  (Quintet  States  of  CO 
Based  on  a  (8e-/3o  +  27i)  SOCI. 


R(a«) 

1  *1* 

2  *1* 

*1“ 

5A 

3.75 

112.664062 

112.619943 

4.0 

112.667000 

112.634461 

112.631212 

4.1 

112.667770 

112.640298 

4.25 

112.668556 

112.646220 

4.5 

112.669242 

112.653411 

4.75 

112.669439 

112.658154 

5.0 

112.669361 

112.661231 

112.660436 

5.25 

112,669134 

112.663194 

5.5 

112.668837 

112.664421 

112.664146 

5.75 

112.668517 

112.665169 

6.0 

112.668202 

112.665609 

112.665594 

112.665P50 

6.25 

112.667908 

112.665855 

112.665915 

6.5 

112.667643 

112.665981 

112.666094 

112.666265 

6.75 

112.667409 

112.666033 

112.666184 

7.0 

112.667206 

112.666042 

112.666221 

112.66633.'< 

7.2 

112.667064 

112.666033 

112.666227 

7.5 

112.666884 

112.666001 

112.666214 

112.666297 

8.0 

112.666653 

112.665934 

112.666168 

112.666228 

9.0 

112.666375 

112.665816 

112.666075 

10.0 

112,666233 

112.665741 

112.666015 

112.666034 

12.0 

112.666115 

112.665669 

112.665959 

20.0 

112.666033 

112.665616 

112.665916 

30.0 

112.666025 

112.665611 

112.665914 

112.665912 

The  con^tations  for  the  and  ^11  states  were  somewhat  more  involved  since 
there  are  two  states  of  each  syrrunetry  which  correspond  asymptotically  to  the  interaction  of 
ground  state  C  and  O  atoms.  Let  us  consider  the  states.  Neither  the  la^  2o2  3<y2  4<y2 
5cfi  Ik^  Irfi  configuration  which  is  dominant  at  large  separations  for  the  1  state  nor 

the  1<^  20^  3<^  4<^  5o  6o  Iti^  2n  dominant  configuration  for  2  ^Z*  correspond 
asymptotically  to  ground  state  carbon  and  oxygen  atoms.  However,  a  proper  average  of 
the  two  mole^ar  configurations  does  correspond  at  long-range  to  ground  state  atoms. 
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Thus,  a  mix  of  one  part  of  the  1  state  and  two  parts  of  the  2  state  dissociates  to 
ground  state  atoms.  We  ran  a  (6e~/2or  +  2ic)  state  averaged  MCSCF  computation  for  the 
^Z'*'  states  widi  the  wei^ting  (1/3)  1  ^Z*^  and  (2/3)  2  At  the  intemuclear  separation 
R  =  30  ao  we  found  the  two  ^Z***  states  to  be  essentially  degenerate  (the  two  states  were 
separated  by  0.028  cm*^)  as  they  should  be.  A  (6er/2a  +  2it)  SOCl  based  on  these  state 
averaged  (SA)  {6er/2a  +  2ic)  MC  output  vectors  yielded  two  roots  which  were  separated 
by  330.29  cm*^.  Obviously,  the  energy  lowering  afforded  by  excitations  to  the  virtual 
space  is  different  (by  about  330  cm*l)  for  the  lowest  two  ^Z*^  states.  We  then  performed  a 
(8e*/3cr  +  2k)  SOCI  based  on  the  (6e~/2a  +  2n)  S  AMC  vectors  and  found  that  the  energy 
sqMUBdan  between  the  two  ^Z***  states  was  reduced  to  90.58  cm*^.  This  SOQ  con^fnitation 

comprised  384,539  CSF;  die  (6e'/2o  +  2k)  SOCI  required  57,898  CSF.  A  SOCI  which 
correlates  ten  dectrons  would  presumably  decrease  die  asymptotic  non-degeneracy  even 
further.  Since  such  calculations  would  have  required  1,M8,672  CSF,  we  did  not  carry 

them  out  The  results  in  TaUe  5  for  the  1 ,2  ^Z'*’  states  correspond  to  a  (8e*/3<J  +  2n)  SOO 
based  on  (6e‘/2a  +  2k)  SAMC  vectors. 

The  situation  for  the  1,2  ^11  states  is  even  more  complicated.  An  equal  mix  of  the 
diree  configurations  1<^  2<^  3<^  4<j2  5o  1«3  2ic2^  1^2  2t^  4<^  5<^  6a  1«2  2n  and 

la^  2e^  3a^  Aa^  5<7  6<t  1r^  2ic  is  needed  in  order  to  assure  asynqxotic  dissociation  to 
ground  state  atoms.  Note  that  the  first  two  configurations  ate  obviously  11  configurations 

while  the  last  one  corresponds  to  the  sectmd  root  of  ^Z'*’  symmetry.  The  (6e*/2a  +  2k) 
SAMC  calculation  at  R3s30 yields  three  roots  which  differ  from  eaA  odier  ^  at  most 
0.231  cm*^  As  with  die  ^Z'*’  states  the  splitting  between  the  two  ^11  states  actually  became 
larger  udien  we  carried  out  a  SOCI  calculatitm.  We  found  that  the  {6trl2a  -f  2r)  SOCI 
required  108,720  CSF,  the  (8e*/3<r  +  2r)  SOCI  required  705,888  CSF,  and  the  (10e*/4o  + 
2ff)  SOCI  would  have  required  3,020,132  CSF  had  we  tried  to  carry  it  out.  We  didn't 
attend  the  latter  congiutation  because  of  the  prohibittve  eiqiense.  In  fact,  we  halted  the 
computations  at  the  (6e*/2<j  +  2r)  SOO  level  based  on  vectors  from  a  50/50  mixture  of 
^Z***  and  SAMC;  these  are  the  results  presented  in  TaUe  6.  Computations  at  the 

(8e'/30  +  2j()  SOO  based  on  the  proper  three-state  SAMC  are  yet  to  be  completed  before 
this  work  is  ready  for  journal  publication. 


10 


TABLE  6.  The  (6e“/2o  +  2k)  SOO  Potential  Energies,  in  Units  of  -Eb,  of  the  1,2 

States  of  CO  Based  on  a  (6e*/2o  +  2ji)  SAMC  Which  is  a  50-50  Mix  of  ^1+ 
and  States. 


R 

E(l«n) 

E(2  sn) 

4.0 

112;575931 

112.546751 

5.0 

112.583429 

112.578824 

5.5 

112.584507 

112.582940 

6.0 

112.585499 

112.584103 

7.0 

112.586348 

112.584423 

8.0 

112.586436 

112.584449 

10.0 

112.586391 

112.584428 

12.0 

112.586372 

112.584438 

20.0 

112.586358 

112.584443 

Di8cmaQn.Qf  Rfiailts  and  Summaty 

Because  of  the  problems  associated  with  the  asymptotic  region,  it  seems  best  to 
discuss  die  calculated  potential  energy  curves  in  terms  of  binding  energy  curves.  We 
obtained  the  binding  energy  corves  listed  in  Table  7  and  shown  in  Figure  2  by  taking  the 

zero  of  energy  figEb  SLUE  to  be  its  molecular  energy  at  Rs30  ao.  The  (6e'/2<7  +  2n) 

results  we  have  for  the  ^11  states  are  only  prelimituuy  since  they  are  not  based  on  the 
correct  SAMC  which  assures  proper  dissociation.  Thus,  we  present  in  Table  7  only  the  1 
state  to  indicate  its  position  relative  to  die  £  and  A  stales.  The  ^11  states  are  {pparendy 

the  most  weakly  bound  of  the  quintets  we  consider,  while  the  1  state  is  the  most 
deeply  bound  of  the  quintets  we  consider.  Based  on  Bussery  and  coworkers  (Ref.  9),  one 
might  have  expected  that  the  state  would  be  slighdy  more  deeply  bound  than  the  2 
state.  However,  our  results  clearly  show  that  2  and  l^A  and  states  are  nearly 

degenerate  over  a  wide  range  of  intemuclear  separations  with  the  1  state  slighdy  less 
strongly  bound.  This  is  rather  remarkable  qualitative  corroboration  of  the  conclusions  one 
can  draw  fitom  Ref.  9. 

Note  diat  the  lowest  ^11  state  is  not  at  all  strongly  bound  in  contrast  to  the  early 

results  of  ONS  (Ref.  13).  While  we  find  the  1  state  to  be  the  most  deeply  bound  of 
the  quintets  we  consider,  it  is  not  nearly  so  deeply  bound  as  was  suggested  by  ONS  (Ref. 
13).  We  have  repeated  dieir  Cl  calculation  for  the  1  ^Z***  state  and  calculated  the  basis  set 
superposition  error  according  to  the  counterpoise  scheme  of  Boys  and  Bemardi  (Ref.  21). 
We  find  their  superposition  error  at  R=4.0  a^to  be  about  2500  cm*l.  This  large  error  is  not 
entirely  unexpected  since,  of  necessity,  ONS  used  a  minimal  basis  set  in  their  pioneering 
work. 
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T  ABLE  7.  Binding  Energies  (in  cm*^)  of  the  Quintet  States  of  CO  Which  Correspond 
Asynqxodcally  to  Ground  State  Caiten  and  Oxygen  Atoms. 


R  (a©) 

1  sx  + 

2  *X  + 

1  5A 

1  sx- 

1  Hi 

3.75 

4.0 

4.1 

4430.8 

-214.0 

-383.0 

+10023.0 

+6836.6 

+5555.6 

+7616.0 

+8692.7 

4.25 

-555.5 

+4255.8 

4.5 

-706.0 

+2677.6 

4.75 

-749.3 

+1636.6 

+642.8 

5.0 

-732.2 

+961.3 

+1202.3 

5.25 

-682.3 

+530.5 

5.5 

-617.2 

+261.2 

+378.8 

+406.3 

5.75 

-546.9 

+97.0 

6.0 

-477.8 

+0.4 

+13.6 

+70.0 

+188.5 

6.25 

-413.3 

-53.6 

-0.4 

6.5 

-355.1 

-81.2 

-77.5 

-39.7 

6.75 

-303.8 

-92.6 

-59.5 

7.0 

-259.2 

-94.6 

-93.5 

-61.6 

+2.2 

7.2 

-228.0 

-92.4 

-68.9 

7.5 

-188.5 

-85.6 

-84.5 

-66.1 

8.0 

-137.8 

-70.9 

-69.4 

-56.2 

-17.1 

9.0 

-76.8 

-45.0 

-35.6 

10.0 

-45.7 

-28.5 

-26.8 

-22.4 

-1.2 

12.0 

-19.5 

-12.7 

-10.1 

-6.5 

20.0 

-1.8 

-1.1 

-0.7 

0.0 

30.0 

0.0 

0.0 

0.0 

0.0 

It  is  clear  diat  the  1  state  of  CO  is  sufficiently  deeply  bound  to  be  thermally 
stable  at  reasonably  low  tenyeratures.  Its  dissociation  energy  of  ~  7S0  cm*^  (about 
1 1(X)  K)  is  about  4  kT  at  0"  C  and  about  14  kT  at  77  K.  The  l^X'*'  state  has  over  ten  bound 
vibrational  energy  levels.  The  2  ^X***  and  1  states,  which  are  neariy  degenerate,  have 
dissociation  energies  under  1(X)  cm*^  (~  140  K);  the  remaining  quintets  are  even  more 
weakly  bound. 

This  suggests  a  way  to  populate  the  1  ^X^  state  selectively.  If  a  mixture  of  groimd 
state  caibon  and  oxygen  atoms  is  kept  at  about  140  K  in  a  magnetic  field  so  that  dieir  spins 
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BINDING  ENERGY  (cm~l) 


aie  kept  aligned,  all  quintets  but  the  1  statr  will  be  thennally  unstable.  Thus,  any  atom 

recombination  will  form  the  1  state  preferentially.  i  i.  "e  may,  of  course,  be  less  subtle 
ways  to  form  the  CO  quintets. 


Figure  2.  Binding  Energy  Curve.s  for  Quintet  States  nf  CQ, 


Much  remains  to  be  learned  about  die  radiative  lifetimes  of  the  quintet  states  and  to 
leam  whether  there  are  ways  to  control  the  various  decay  processes,  before  we  can  assess 
folly  the  potential  of  high-spin  CO  as  an  early  storage  mole^e. 
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APPENDIX 


Semi-Empirical  Deteimination  of  the  Long-Range  Molecular  States 
Dissociating  to  C(3Pj)  +  0(3Pj) 


The  following  material  is  an  incon^lete  semi-final  draft  of  a  manuscript  we  have 
prepared  to  submit  for  publication  in  Chemical  Physics.  It  is  intended  only  to  provide  the 
det^  necessary  to  tqtptedate  die  contents  of  the  fim  part  of  die  technical  discussion  in  this 
rqiort. 
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I.  Introduction 


The  thetuedcal  or  senu-eii^>iiical  treatment  of  the  Itmg-range  interaction  of  a  pair  of  atoms  is  best 
developed  for  systems  in  which  the  interacting  fragments  are  closed-shell  [1]  or  else  where  one  or  both 
fragments  differ  from  a  closed  shell  electronic  configuration  by  a  single  electron.  For  example,  the  semi- 
empirical  model  of  Busseiy  and  Aubert-Frecon  (BAIO  has  been  quite  successfrd  in  either  piloting  or 
reproducing  the  very  long  range  interaction  energy  curves  for  alkali  or  halogen  diatomic  molecules  [2-4]  as 
judged  by  comparison  with  experimental  spectroscopic  or  exacting  theoretic  ab  initio  results.  The 
situaricm  is  not  nearly  so  sanguine  for  the  interactitm  of  mote  general  open-shell  fragments,  however.  It  is 
our  objective  in  tfiis  piper  to  e3q>lote  the  extension  of  the  BAF  nnodel  to  the  interaction  of  open-shell 
fragments. 

We  choose  presently  to  study  the  interaction  of  ground-state  open-shell  C  and  O  atoms.  Our 
motivation  tot  this  dioioe  stems  in  part  from  the  interest  in  characterizing  the  weakly  bound  quintet 
mdecular  states  resulting  in  the  interacticm  of  atomic  C  (^P)  +  0(3P)  as  possible  eitergy  storage  states. 
Furthermore,  there  are  new  ab  initio  cmipotatimis  [S]  and  spectr^opic  observations  [6]  in  progress  with 
whidi  to  conpare  our  results.  Toward  this  end,  w»  present  here  die  long-range  interaction  coefficients  for 
the  dectrostatic  and  dispersion  energies  and  die  resulting  potential  energy  curves  which  correspond  to  dw 
C  (3P)  +  0(3P)  interacrion. 

hi  Section  n  we  review  the  BAF  mediod  as  tpplied  to  the  CO  molecule.  We  consider  both  the 

[A,S]  and  [JJ]  angular  momentum  coiqding  schemes  which  are  appropriate,  respectively,  to  the  cases 
where  the  atomic  spin-orbit  energy  is  small  or  larro  conpared  to  the  interaction  energy.  We  first  build  ip 
the  symmetry-adipted  tndecules  wavefiinctions  for  atrans  in  dieir  ground  term  or  ground  level  for  these 
two  coipling  cases.  First  order  perturbation  theory  enables  os  to  compote  the  quadnipole-<{uadnpole 
energy  terms  leading  to  the  C3  coefficients.  Second  order  perturbation  theory  gives  rise  to  die  induction 
and  dispersion  terms  diarac^ued  by  the  €5,  Cg,...  coefficients. 

In  Sectkm  in  we  discuss  details  of  the  calculations.  There  we  demonstrate  die  sensitivity  of  our 
results  to  iiput  parameters  of  die  modd  such  as  the  value  of  <fi>,  the  mean  square  radid  position  of  the 
electrons  in  the  ground  stale  atoms  and  the  line  strengths  of  electric  d^le  transitions  of  the  atoms. 

In  Sections  IV  and  V  we  examine  both  HuikTs  case  "a"  and  case  "c”  interaction  energy  curves  for 

intemuclear  separations  ^  j  suggested  by  the  Le  Roy  criterion  [6]: 

(for  CX)  R  ~  6.7  bdir).  In  the  (A,S)  coupling  case  the  present  model  suggests  somewhat  weaker  binding 
than  obtained  in  the  ob  initio  conputations  of  Rosenknuitz  and  coworkets  [5].  Doe  to  the  strong  spin-orbit 
irtteraction  in  oxygen,  however,  we  conclude  that  the  (JJ)  coipling  scheme  most  be  used  for  an  adequate 
descri^on  of  di^  very  long-range  interactions. 


n.  Modd 

The  mediod  we  use  in  this  work  is  basically  die  same  as  that  already  described  and  used  for  die 
hdogen  and  aHcali  dimers  [2-4].  The  present  summary  is  conseqoendy  brief.  We  investigate  the 
molmnilar  stales  due  to  the  interaction  of  the  ground  state  atoms  C  and  O  for  intemuclear  distances 
sufficiently  large  so  that  the  overlap  between  the  atomic  wavefonctions  is  presumed  to  be  negligible. 
Conseqoendy,  we  neglect  ovethp  and  exchange  effects.  In  the  usud  perturbation  theory  treatment  up  to 
second  order  we  invd»  the  multipolar  expansion  of  the  perturbation  operator: 
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where 


V  = 


ial  jsl  R 


i+J+1 


Vij(a,b)  =  (-)i  ^  dm(i,j)Q;„(a)Qi^(b) 


and 


dm(i.  j)  =  (i  +  j)  t  ((i-m) !  (i  +  m)  I  (j  -  m) !  (j  +  m) 


!)■ 


1 

2 


(r,  9,  ^  are  the  spherical  polar  coordinates). 


The  spin-orbit  separation  of  43.4  cm'l  for  C  and  227  cm*^  for  O  validates  the  LS  coupling  for  the 
atcmis  and  the  Hand's  case  co(q)ljng  ''c”  for  die  molecule.  The  molecular  stales  may  then  be  classified 

following  the  rqnesentation  ixdiere  Q  *  t  Mj^  Mr)  I  is  the  absolute  value  of  the  projection  of  the  total 
electronic  angular  momentum  along  the  intemuclearaods.  The  superscripts  ±  are  related  to  the  states 
n  s  0.  In  this  representation  the  interaction  between  C(^Pj)  -f  O(^Pj)  gives  rise  to  49  molecular  states; 


states 

0+ 

o-.i 

o-.i 

0+  0+,  0-.  1,  1,  2 
0+.  1.  2 
0+  1.  2 

0+  0-,  0-.  1,  1, 1,  2,  2,  3 
0+  0-,  0-,  1. 1,  1,  2,  2,  3 
0+  0+,  0+,  0-,  0-,  1,  1,  1,  2,  2,  2,  3,  3,  4 


dissociation 
C(3Po)  +  0(3Po) 
C(3Po)  +  0(3Pi) 
C(3p,)  +  O(3po) 
C(3p,)  +  CX3Pi) 
C(3Po)  +  0(3P2) 
C(3P2)  +0(3Po) 
C(3Pi)  +  CX3P2) 
C(3P2)  +  0(3Pi> 
C(3P2)  +  0(3P2) 


The  zero  order  molecular  functions  needed  in  the  perturbation  procedure  are  built  up  as  simple  products  of 
non-overii^iping  atomic  wavefunctions  labeled  a  and  b.  Following  Chang's  [7]  prcMredure  each  function 
may  be  written  as: 


Ov 

Vo=  '  JkMfc>b 

ksl 


The  basis  fimctions  yu  are  presented  in  Appendix  I  for  the  49  molecular  states. 
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We  also  investigate  the  interaction  between  the  two  ground  state  atoms  C  and  O  where  we  neglect 
'  the  ^in-orbit  effect  wMiin  each  atom.  The  corresponding  molecular  states  may  then  be  classified 

following  the  2Sfl^±  representation  where  S  is  the  total  spin,  A  is  the  projection  of  the  total  orbital  angular 
momentum  on  the  intemuclear  axis,  the  superscr^  ±  denoting  the  symmetry  pr<^)erties  of  the  coordinate 
wave  functions.  In  this  representation  (Hund's  case  coupling  'V)  the  interaction  gives  rise  to  18 

rntdecular  states:  two  stales  of  each  symmetry  two  states  of  each  symmetry  one  state  of 

each  symmetry  one  state  of  each  symmetry  dissociate  to  C(3p)  +  0(3p).  The  molecular 

basis  functions  dnived  in  Chang's  [7]  procedure  may  be  written  as: 

no 

yo=  I  Sk  Lk  Msk  MLk  >a  •  Ml^  >b 

ksi 

The  e:q>anded  f<»ms  of  the  yo  ate  presented  in  Ai^rendix  n  for  the  18  tiwlecular  states. 

The  terms  in  the  matrix  elemertts  of  the  interaction  energy  which  depend  only  on  the  angular  part  of 
the  atomic  wave  funcdons  (quantities  A  and  C  in  equations  (7, 15, 25, 30)  of  ref  [4]  remain  valid  for  CO 
while  die  terms  which  dep^  on  the  radial  part  of  dw  atomic  wavefutictions  take  on  specific  forms.  In  the 
paiticolar  case  of  the  carbon-oxygen  interaction  in  their  ground  term  ^P,  the  initial  cor^guration  to  be 
consideied  is  Is^  2s^  2^  for  the  carbon  atom  and  Is^  2^  2^  for  the  oxygen  atom.  For  bodi  of  them  the 
optical  electron  is  included  in  one  group  of  equivalent  electrons  p2  or  ^  and  has  to  be  decoupled  from  the 
others.  In  the  framework  of  the  fractional  parentage  scheme,  the  resulting  expiessions  for  tte  various 
radial  quantities  ate  die  following: 

a)  for  die  electrostatic  energy  terms: 

Ri»<(l“)L.IIQi|l(l»)L.> 

L,  V  V  }  (o  0  J)<nllr‘lnl>. 

is  given  by  a  relation  similar  to  Eq(l),  where  primed  atomic  quantum  numbers  ate  used. 

SL 

The  (piantities  are  the  fractional  parentage  coefBcietits  for  the  |  ( l'‘‘l[SiLi]l)l°SL} 

configuration  [8].  The  sum  is  over  the  various  spectral  terms  SiLi  of  the  configuration  Prom  the 
selection  rules  dw  to  the  non-zero  conditions  of  the  3j-Wigner  coefficients,  the  only  allowed  values  are 
t=i's  2  which  give  rise  to  a  C5  contribution  in  the  electrostatic  energy  term: 

Ee*««  a  E22(t)  A)  *  E5  (u,A.)/R5 
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b)  for  the  dispersion  energy  terms: 


R 


r  ij'  ^  S 


<Y.  L,  II  Q‘  II  n  L'iXri  L"  II  Qi'  II  r,  Li> 


(2) 


»  <K  Li,  II  Qj  II  II  Qi'  II  It  Lt>  [aE  (n"  1")  +  AE  (n't  I't)] 

where  Y  includes  all  the  quantum  numbers  necessary  to  represent  die  electronic  state.  The  sum  is  carried 
over  all  the  electronic  excited  states  of  each  atom.  Jh  the  particular  case  of  the  carbon-oxygen  interaction 
duee  types  of  excitatitm  can  occur 

i).  Excitation  1“ - l“-l  (n’T") 

The  excited  optical  electron  is  one  of  the  (1<>)  electron  group.  Only  one  group  of  equivalent 
electrons  needs  to  be  (XMisidered.  The  matrix  elements  can  be  written: 

<(in)  La  II  Qi  II  l“-i[SiLi]  n''V;  L'><in-l[SiLi]  n';  T'  L';  II  Qi'  II  (!»)  L,> 
'"[^stLir  (2U+1)(2L;'+1)(21+1)(21;'+1)  (2a) 


1 


l1  1*  }(  0  0  0*  )(  0*  0  0  I  f*  I  n  ,  1'^  <n',  I''  I  ri  I  nl> 
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ii)  Excitation  1“  1'*’ - 1*”"^ 

The  excitatitm  occurs  between  two  groups  of  equivalent  electrons.  The  matrix  elements  are  given 
by: 

<in[SiLi]  I'P  [S2L2]  SaLa  H  Q'H  [Sf  L{]  1^]  L'pxS^  L'^  II  Qi'  II  Sa  La> 


(2^+l)(-faLa 


f  u  r  L3  If  4  r  L3 

(2L3  +  1)  (2L3  +  1)  (2La  +  1)  (2L  a  +  L)  (21i  +  1)(21  +  L)  (21'+  2)  |  l;'  j{  L2  L';  4 

J  Li  La  L2  If  L3  U'  L2  U  1  L2  Lf  1  f  r  4  41/1  i  1'  \  /  1'  i'  1  \ 

1  L'a'  L3  i  n  La  Li  i'  A  L3  L  i  M  Li  1  r  jU  0  0  M  0  0  0  ^ 


(2b) 


<nl  I  ri  I  n'lXn'l"!  ri'  I  nl> 
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iu)  Excitation  1“  1^- - in-l  n*'l" 

Two  groups  of  equivalent  electrons  have  to  be  considered  but  only  one  is  modified.  For  that 
particular  case,  we  give  tire  general  expressions  for  the  multi^lar  transition  moments. 


<(1“)  [SiLi]  I'*’  [S2L2]  SaLa  II  Qi  II  S3L3  n'^  V;  L';>  (2c) 

X  <S3L3  n'.'  V;  L'.ll  Qi'  II  1«'  [S2L1]  IT  [S2L2]  SaLa> 


(2S3  +  1)(2Si  +  1)^ 


,0U*1)(2L'.-+1){4  {s  ■'  J,) 


'"[“swT 


S3 
Si  S2 


X  <nl  I  r*  I  n-  l'i>)<n''  T'  I  r‘  I  nl> 


F(V  the  contribution  of  atom  b,  expressitms  similar  to  Eqs.  2(a-c)  may  be  derived  using  the  substitutions 
L''  — >  L'^ ,  La  — >  Lb.  The  selection  roles  due  to  the  non-zero  conditions  of  the  3j-Wigner  coefficients  lead 
to: 

00 

E<fisp  (\),X)  «  X  E2n  (V.  X)  /  R2  “ 
n»3 

where  2n  =  i  +  i"+  j  +  j".  Thus,  only  even  powers  of  R  are  involved  in  die  multipolar  expansion  of  the 
dispersion  energy. 


ni.  Calculations. 

For  both  the  electrostatic  and  disposion  energy  terms  whether  or  not  we  include  ^in-oibit  effects, 
we  obtain  each  matrix  element  as  a  sum  of  products  of  purely  angular  terms  (eqs.  7,  IS,  2S,  27  of  ref.  [4]) 
with  radial  ones  (eq.  S.  1  and  2).  Angular  terms  are  conqnited  ex^y  firom  die  Wigner  3j  and  6j 
coefficients.  The  numerical  evaluation  of  the  radial  terms  requires  die  kiKiwledge  of  the  r^al 
wavefimctions  for  the  initial  state  as  well  as  for  all  the  excited  states. 

For  the  investigation  of  the  carbon-oxygen  interaction,  we  have  chosen  the  following  (qiproximate 
wavefiitictions: 

a)  The  Carbon  and  Oxygen  2p  and  2s  orbital  wave  functions  displayed  in  i^ipendix  IU  have 
been  determined  from  SCF  calculations  using  the  ALCHEMY  program  [9].  In  order  to  dwck  die  quality 
of  die  present  ground  state  wave  functions,  we  compare  in  Table  1  present  values  with  other  ab  initio 
calculations  [10-14]  of  <t2>.  Good  agreement  is  found  with  Froese-Fischer's  Hartree-Fock  calculations 
[12].  But  our  values  ate  small  conqiared  to  Desclaux'  Dirac-Fock  calculation  [14]  or  other  recently 
published  data  [1 1]  for  carbon  as  well  as  for  oxygen  atoms  and  result  in  a  smaller  value  of  R.  To  improve 
the  C5  evaluation,  the  Desclaux  results  for  <r2>  have  been  used  in  our  conq)utations. 

b)  The  excited  wavefvmctions  have  been  determined  in  the  Coulomb  qjproximatioa  They  are 

described  by  a  Whittaker  function  where  the  effective  quantum  number  u  and  the  effective  charge  z*  are 
determined  in  the  hydrogenic  model.  That  is. 
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_  _  i  _  i 

\)  =  l/[2(E[SiLi]  -  E[SiLi]  nlSL)pandz*  =  n[2(E[SiLi]  -  E[SiLil  nlSL)]2. 

Here,  E[SiLi]nlSL  is  the  aven^  value  of  the  experimental  electronic  state  energy  defined  by; 

EiSiLi]nlSL=^  (2J+1)  E[SiLi]nlSU  (2J+1) 

J  J 

and  E[SiLi]  is  the  average  energy  value  of  die  2Si  +  ILi  terms  which  is  equal  to  or  larger  than  die 
ionization  potential  of  the  neutral  atom. 

For  die  n"s  series  the  quantum  defect  p=n-t)  is  roughly  constant,  so  we  represent  the  excited  state 
by  a  Whittaker  function  with  an  effective  quantum  number.  For  the  n^d  series  the  effective  charge  z* 
keeps  roughly  the  same  value  so  we  repiesoit  the  excited  stales  by  a  Whittaker  function  with  an  effective 

charge.  For  highly  excited  states  not  known  experimentally,  we  use  a  Whittaker  function  where  V  =  n  and 
z*  s  1  which  makn  it  equivalent  to  a  hydrogenic  function.  We  have  limited  our  evaluation  to  the 
deteimination  of  die  Q  coefficient  for  i^di  only  the  n''s  and  n"d  series  are necessary. 

In  TaUe  2a  and  2b  we  compare  the  line  strengdis  for  some  carbon  and  oxygen  electric  dipole 
transitions  computed  in  this  work  (using  the  apjxoximations  ^iriiich  we  have  just  d^cribed)  to  odier 
thetxetical  and  experimental  values.  The  line  strengths  (in  a.u.)  are  related  to  oscdlator  strengths  or 
emission  transitkm  probabilities  by  the  relations: 

g.Xfik/303.8  and  S  *  g^X3Au/2.0264  x  lO^ 


where  g.  and  ^  ate  the  statistical  weight  of  the  lower  and  iqiper  state  respectively,  X.  is  the  transition 

wavelength  in  (in  A),  fjk  is  the  absorption  oscillator  strength  and  Au  is  the  emission  transition  probability. 
Furthennote,  we  have: 

S  =  (2s+l)  (<Y  L  IQl  |l  Y’L"»2 

where  s  is  the  atomic  spin  quantum  number  and  <Y  L  P  P  Y'L">  is  the  dipole  transition  moment  given  by 
Eqs.  2a,  2b  and  2c  using  i=i'=l. 

Table  2  lists  the  transition  moments  we  have  calculated  together  with  values  taken  finm  the 
literature  (15  -  41).  We  have  found  few  literature  values  for  the  2s^  ^  -  2s  2p^"l "  transitions  in  carbon 
wM  which  to  ootttptae  our  calculated  results.  No  transition  moment  or  related  quantity  determination  has 
been  reported  very  recently  for  diat  atortt  Overall,  our  tqtptoximate  values  agree  best  with  the 
e;q)erimental  results  of  Boldt  [19].  The  agreement  is  go^  for  the  [2po]  3d  excitation  while  our 
calculated  values  ate  snail  congtated  to  Boldt's  fcv  the  [2po]  3s  and  [2po]  3d  transitions.  It  should 
be  noted  in  Table  2-a  that  diere  often  exist  serious  dtscrepatKies  between  the  various  experimental  values 
of  the  line  strerqpdis,  so  disagreement  with  experiment  is  not  necessarily  an  indication  of  inadequacies  in 
our  calculations. 

We  have  taken  for  carbon  the  most  recent  and  presunufoly  most  accurate  theoretical  transition 
nKxnents  given  by  Froese-Fischer  and  Saha  [38]  who  used  the  MCSCT  method  with  Breit-Pauli 
corrections.  A  cotrqtarison  of  their  transition  moment  values  (listed  in  the  column  labeled  FF)  in  Table  2-a 
with  other  calculated  values  shows  that  cortelation  effects  ate  sf^Nuentiy  quite  ingxntant  for  the  transition 
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between  2S^  2p^  and  2S  2P^  ^S,  ^P  and  ^D.  Neverdieless,  discrepancies  between  dieoretical  and 
experimental  results  still  persist  and  are  in  some  cases  important,  especially  for  the  2s-2p3  3p  transition.  In 
or^r  to  see  the  influence  of  this  choice,  we  have  also  ccnnputed  the  C5  coefficients  using  the  results  of 
Nicolaides'  many  body  calculations  [34]  for  the  2s  ^  transition  moments  which  are  seen  to  be  in  good 
agreement  with  expeiimaital  values. 

In  Table  2-b,  we  compare  our  present  values  for  the  electric  dipole  moments  for  the  2s  2p^  nT  and 
2s  2p^  transitions  in  oxygen  widi  other  theoretical  and  erqrerimental  dW  Perhaps  the  most  reli^le 
confuted  values  are  tfrase  of  Froese-Fischer  [2S]  who  tilled  a  MCSCF  method  widi  Breit-Pauli 
correcdmis  for  the  relativistic  effects.  Our  values  are  too  small  for  die  term  while  Froese-Fischer's 
value  for  the  3s  term  is  dose  to  the  experimental  one.  Our  present  results  are  dose  to  those  due  to 
Armstrong  [26]  whose  extended  Coulomb  approximation  calculations  are  similar  to  ours.  For  the  [^D^  3s 
3d  and  [2P^  38"  3p  terms,  our  values  are  too  large  while  Froese-Fischer’s  compare  well  with  experiments 
for  die  3d  teim.  All  die  theoretical  values  for  die  3p  term  are  very  low  crmqiared  to  dw  single  experimental 
value.  There  may  be  a  large  error  in  eiqierimental  results  for  the  ^P^  3s"3p  transition  due  to  the 
autoionizing  chanuxer  of  foat  term.  Alro,  Pradhan  and  Sanph  [27]  have  computed  the  oscillator  strengths 
for  diese  transitions  based  on  a  close  coupling  mediod  in  the  frozen-core  qiproximation.  Their  values  are 
larger  than  those  of  Froese-Fischer.  For  the  n”s”  and  n”d’’  series,  our  values  are  similar  to  Pradhan  and 
Sarsph's  [27]  and  Chung's  [28]  results  for  the  3$  and  3d  terms.  Chung's  calculations  are  based  on  the 
Hattree-Fock  single-configuration  wavefunctions  with  a  Is^  2s^  frozen  core  (our  core  is  not  frozen).  Our 
values  are  between  those  of  the  two  authors  for  the  n"s"  series  and  very  close  to  (Thung's  values  for  the 
n"d"  series,  but  lower  than  die  experimental  data  [18]  for  bodi  die  n"s"  and  n"d"  series.  Our  values 
decrease  widi  n"  more  quickly  dian  do  die  experimental  ones.  This  fact  may  have  a  non-negligiUe 
influence  on  the  C^  evaluation.  A  large  discrqiaicy  can  be  seen  for  die  (^I^]3d"  3p  Une  strengths 
between  experimental  and  themetkal  determinations.  The  Froese-Fisdier  value  of  the  transition 
probability  for  the  2s  2p3  3p  transition  in  oxygen  and  is  found  to  lie  between  experimental  and  odier 
diemetical  results.  This  value  has  been  iiKlu^  in  our  computation  as  well  as  that  of  Nicolrides  [34] 
which  gives  rise  to  different  values.  From  these  ctmqiarisoiis  results  the  fact  that  for  the  first  electric 
dipole  transition  moments  (at  least  the  3s,  3d  excitations),  our  qiproximate  scheme  does  not  give  reliable 
results.  The  correlation  effKXs  are  preponderaiit  in  these  carbon  and  oxygen  transitions  and  more  conqilex 
calculations  dian  the  present  ones  are  necessary. 

The  evaluation  of  die  quantities  requires  furdiermore  the  knowledge  of  all  the  excited  state 

energies  which  we  take  as  the  experimental  ones  [35].  We  have  included  20  excited  states  in  the 
summation  in  Eq(2).  Higher  excited  states  do  not  make  significant  contributions  to  the  long  range 
coefBcient  value. 


IV.  Results  and  Discussitm 

We  have  performed  the  perturbation  theory  calculation  of  the  first  order  quadrupole-quadtupole 
interactions  as  well  as  second  o^r  induced  dqxile-induced  dqxile  interactions  for  the  system 
C(3Pj)  -h  0(^j).  We  write  each  interaction  energy  matrix  element  as: 

Ete(1>.X)  =  8ttXiEdh(«) 


where  Edis  (h)  is  the  dissociation  energy  for  the  molecular  state  v.  Let  us  first  consider  the  case  where  we 
neglect  spin-orbit  effects. 
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Dae  to  die  neglect  of  the  exchange  contributim  the  various  E5  and  E<$  terms  are  identical  for  the 
mult^ets  1, 3, 5  of  each  symmetiy.  Thus,  the  18  x  18  interaction  energy  matrix  reduces  to  4 

submatrices.  There  is  one  of  order  2  for  each  of  symmetries  h*  and  11  and  one  of  "order  1"  for  each  of 

the  states  A  and  (See  ^rpendix  IV  for  explicit  expressions  of  the  £5  and  £5  matrices).  For  the  states 

Z*  and  n,  the  corresponding  submatrices  have  to  be  diagonalized  for  various  values  of  die  intemuclear 

sqparation  R,  thus  providing  an  interaction  energy  curve  £u  (R)  f(»  each  state  M .  These  carves  are  shown 
in  Fig.  1. 

We  have  evahiated  die  energy  eiqnessions  for  values  of  R  starting  from  7  a.a.  and  increasing  by 
steps  of  0. 1  a.a.  The  lower  limit  is  suggested  by  the  Le  Roy  criterion  R  >  2^<r^>  +  <r^p.  By  a 
least  squares  fit  of  these  energy  values  for  each  state  to  the  eiqnession: 


E^(R) 


Cs(t))  .  C6(t)) 

TT*  r6 


+  £di8(>)) 


allows  us  to  determine  fitted  values  for  die  two  long  range  coefficients  C5  and  €5  for  each  molecular  state 
£'''andrL  These  values  are  presented  in  TaUe  3. 

hi  the  fits  we  have  induded  increasing  values  of  R  until  die  dissociation  limit  was  reached  (at  a 
precision  of  O.QS  cm*^).  Note  that  under  the  present  assungition  of  neglecting  spin-orbit  effects  within 
each  atom,  the  eighteen  states  carespond  to  a  unique  averaged  dissociation  l^t.  The  values  of  the  Q 

coefficients  for  the  Z*,  £*,  11  and  A  molecular  states  are  all  within  a  factor  of  two  of  each  other.  Thus,  the 
main  difference  between  dK  various  interaction  energies  results  fiom  their  different  Cs  coefficients. 

The  different  values  which  are  listed  in  Tdde  3  show  the  sensitivity  of  these  quantities  to  the 
choice  of  the  2s^  ^  -  2s  transition  moments  and  the  choice  of  the  first  3s  or  3d  atomic  transition 
moments.  We  have  compi^  die  €5  [#1]  estimations  using  the  qiproximations  discussed  in  Sectioil  in 
for  all  di^le  transition  moments  and  Nicolakfes'  estimates  for  the  2s^  2p°  -  2s  transition  moments 
for  carbon  ^S,  ^P,  terms  and  the  <^gen  term  [34].  We  have  congnited  the  Q  [^1  estimations 

using  these  same  iqiproxiniations  as  in  #1  except  that  here  we  use  Froese-Fischer's  estimations  for 
2s^  ^  -  2s  2^1  transition  maments  in  carbon  [38]  and  oxygen  [25].  We  have  conqxited  the  Q  [#3] 
estimations  using  values  for  the  first  dqiole  transition  moments  (3s  and  3d  excitations)  obtiuned  fi<^ 
e:qietimetu,  correlated  connotations  (selected  fiom  Table  2a  and  2b)  or  our  own  estimates  and  Froese- 
Fischer’s  results  for  die  2s^  ^  -  2s  transitions  nuHnents  togedier  with  the  q)p>^<”(intations  used  in 
#1  for  higher  excited  states.  It  is  evid^  that  die  coefficients  are  very  sensitive  to  die  determination  of 
the  first  dipole  transition  moments  wiudi  contribute  60%  of  the  total  value  (note  that  this  contribution  is 
about  90%  for  alkali  dimers). 

A  Q  estimation  iMsed  on  Unsold's  rqiproximation  has  been  given  by  Qiang  [7].  He  estimates 
=  20.8  a.u.  and  C5  » 10.65  a.a.  for  the  states.  Our  C5  value  for  these  molecalar  states  is  in  good 
agreement  with  Oiang's  determination  while  our  value  is  much  lower  than  his  estimation.  The 

2 

accuracy  of  the  C5  depends  mainly  on  die  accuracy  of  the  values  used  for  the  averaged  quantities  <r^>, 

2 

<rQ>  which  apparendy  is  good  (see  TaUe  1).  For  the  coefficients,  the  accuracy  of  the  present 
determination  is  limited  by  the  accuracy  of  the  values  of  the  transiticm  moments  we  use.  Note  that  our 
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calculated  values  of  these  qpiantides  were  seen  to  be  in  good  agreement  with  (but  generally  smaller  than) 
those  available  fiom  experiment  or  other  theorerical  calculations.  We  note  also  that  experimental  values, 
when  they  are  available,  are  sometimes  given  with  rather  large  uncertainties.  There  is  a  particularly 
significant  q)read  in  the  values  for  die  carbon  transition  2s  2p^  ,  for  example;  di^rent  values 

S  =  12.47  [39]  or  S  =  3.62  [34]  are  at  our  disposal.  The  Coulomb  £q>proximation  gives  reasonable 
transition  moments  for  Carbon  and  Oxygen  excited  states  as  judged  by  comparison  with  the  available 
eiqierimental  values,  .^jpaiendy,  only  the  first  few  transition  moments  are  sensibly  af^ed  by  correlation 
effects.  Conse^endy,  obtaining  improved  C$  values  depends  sensibly  on  obtaining  more  reliable  values 
of  die  first  few  transition  moments  from  either  dieory  or  experiment.  We  have  found  dw  contribution  of 
the  di^iersion  energy  to  the  long-range  interaction  oiergy  to  be  relatively  small  (C^  less  than  20  a.u.) 

crenpaied  to  the  quadrupole-quadiupole  interaction  energy  for  the  1  and  2 11  states.  Indeed,  we  have  at 
R  =  7  a.u.  (vdiich  is  die  lower  limit  of  the  "pure"  long-range  interactions)  C6/(R  C5)  =  0.1 1  for  the 

1  1>3,S£+  states  and  0.16  for  the  2  states.  For  the  states,  the  dispersion  energy  due  to  the 
induced  dipole-induced  d^ole  interaction  is  of  the  same  order  as  the  electrostatic  one.  We  have 

C^C5*R)  —  ~  0.65.  So,  for  the  1*3,5^  states,  a  very  accurate  C6  coefficient  would  be  more  significant 
than  it  is  to  other  states. 

No  eiqierimental  electronic  energy  curves  of  the  CO  molecule  are  known  up  to  now  for  intemuclear 
separadtxis  as  large  as  7  a.u.  Recendy,  however,  Rosenkrantz,  Bohr  and  Konowalow  (RBK)  [S]  have 
peifonned  multi-refeienoe  second  order  configuration  interaction  calculations  on  the  quintet  states  of  CO 
which  c<XTespond  asyngitotically  to  the  intoacdon  of  ground  state  C  and  O  atoms.  With  the  exception  of 

the  2  3n  state  which  is  essentially  repulsive  each  of  the  quintets  is  relatively  weakly  bound  with  a  potential 
minimum  which  lies  at  an  intemuclear  separation  R  >  4.75  ao  (By  conqiarison,  the  strongly  bound  ground 
state  has  its  potential  minimum  at  about  2  ao).  RBK  find  somewhat  greater  bnuling  than  we  do.  If  we 
assume  diat  die  RBK  curves  are  neatly  exact,  it  is  evidtmt  that  the  ovetiap  and  exchange  effects  have  a 
non-negligible  contribution  to  the  interaction  energy  well  beyond  the  limit  (6.7  a.u.)  given  by  the  Le  Roy 
criterion.  The  difference  between  the  "pure"  van  ^r  Waals  interactions  isiiich  we  calculate  here  and  the 
conqxited  energy  curves  of  RBK  give  us  an  estimate  of  exchange  and  danqnng  e^cts  and  the  magnitude 
of  the  dispersion  terms  we  have  neglected  to  eadi  molecular  state.  In  fact,  at  such  intemuclear 
sqiarations,  the  interaction  energy  to  all  die  electronic  states  dissociating  to  C(3p)  +  0(3p)  (except, 

possibly,  for  the  Z*  states)  is  small  compared  with  the  spin-orbit  raergy  of  the  oxygen  3p  term.  So,  the 
spin-orbit  interaction  must  be  taken  into  account.  As  a  fim  {pproximation,  we  include  the  spin-orbit 
energy  on  each  atom  and  present  the  electronic  eneigy  curves  in  the  so-called  Hund's  case  "c"  coupling 
scheme  which  gives  rise  to  49  molecular  states. 

The  49  X  49  energy  matrix  Eim  is  factored  into  6  submatrices  of  order  10  to  the  symmetry  0+  of 
order  9  for  the  symmetry  O',  of  order  15  for  the  symmetry  1,  of  order  10  for  the  symmetry  2,  of  order  4 
fordiesymmetty  3,  andof  "order  1"  tothestate4.  The  Es  and  Ee  submatrices  are  presented  in  Ai^ndix 

V  (a-e)  for  the  O**,  0, 1, 2  and  3  symmetries  respectively.  For  the  4  symmetry  molecular  state,  and  oidy 
for  this  symmetry,  we  have  C3  =  E5,  €5  =  E$  and  €5(4)  =  C5  (1*3,5^)^  C6(4)  -  Cg  (13.5^),  por  each 
symmetry  we  diagonalize  the  corresponding  submatrix  for  numerous  values  of  R  providing  Eo(R)  for  each 
state  v.  Graphic  representation  of  the  49  molecular  states  are  shown  in  Fig.  2  (a-e). 

Some  long-range  extrema  are  seen  for  the  O^,  O',  1, 2  and  3  symmetry  states.  We  present  in 
Table  4  their  position  (in  a.u.)  and  interaction  energy  (in  cm'^).  The  nindi  state  of  symmetry  "1"  and  die 
eighth  state  of  symmetry  "2"  have  a  well  depth  of  1 1.8  and  15.5  cm*^  respectively.  The  other  states  which 
exhibit  a  minimum  have  a  wdl  depth  less  than  5  cm*^;  this  is  smaller  than  the  precision  we  have  on  each 
molecular  state  and  is,  consequently,  not  significant.  At  shorter  intemuclear  separations  the  (J  J)  coupling 
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molecalar  states  give  rise  to  the  (A,S)  coupling  states.  But,  as  is  suggested  by  the  present  results  shown 
in  Fig.  2,  the  (A,S)  coiq>ling  scheme  is  not  accurate  enough  even  at  R  =  7  a.u.  to  represent  the  long-range 

interactions  of  C(3Pj)  +  0{3Pj)  except,  possibly,  for  those  leading  to  the  2+  states.  We  give  in  Fig.  3  a 
correlation  diagram  between  Hund's  case  "c"  and  Hund's  case  "a"  molecular  states.  By  comparing  the 
curves  for  the  various  states  described  by  die  quantum  number  A  or  £1,  it  appears  that  the  "4"  symmetry 
state  comes  from  the  state.  We  note  more  paiticulariy  that  die  first  "3"  symmetry  molecular  state  comes 
from  the  Estate,  die  second  comes  from  the  ^A,  the  third  comes  fiom  the  1  ^11  state  and  the  fourth  comes 
from  the  2  state. 

V.  Summary  and  Conclusion 

The  present  paper  reports  estimations  of  potential  energy  curves  and  long-range  coefficients  of  the 
mcdecular  states  of  CO  dissociating  to  C(^Pj)  +  0(^Pj).  These  interactions  give  rise  to  singlet,  triplet  and 
quintet  1,*,  £*,  11  and  A  molecular  states  in  (A,S)  coupling  or  0+,  ()-,  1, 2, 3  and  4  symmetry  states  in 
(J,J)  coiqding.  They  result  in  a  quadrqxile-quadrupole  interaction  (C5/R^)  energy  term  due  to  the 
permanent  quadnqwle  of  each  carbon  and  oxygen  atom  in  their  ground  state  configuration,  plus  dispersion 
interactions  (C6/R^  +  Cg/R^  +....)  energy  terms.  Since  the  requisite  input  data  is  so  uncertain  we  corrqxite 
ordy  the  induced  di^le-induced  dqxile  energy  terms  (Cg/R^)  for  each  state.  The  C6R‘^  interactions  are 
quasi  independent  of  the  molecnlar  state  since  C^  is  nearly  the  same  for  die  1,*,  2',  11  and  A  states. 
Futthermore,  the  C^/R^  mergy  contribution  is  small  relative  to  die  quadnqxde-quadruprde  energy  for  the 
1 2'*’  and  2 11  states  but  is  of  the  same  order  for  die  1 A  states  while  the  electrostatic  energy  vanishes  (C5 

s  0)  for  die  2  2'*',  1 2*  and  1 11  molecular  states.  The  C5  coefficients  are  given  with  a  good  accuracy 
while  die  accuracy  the  C6  coefficients  dqiends  mainly  on  the  uncertain  quality  of  the  electric  dqxile 
moments,  e^iecii^y  those  evaluated  in  the  Coulomb  approximation.  The  evaluation  of  the  dispersion 
energy  is  die  most  significant  shortcamnig  of  the  BAF  tqiproach  as  ap^ied  to  long-range  CO  interactions. 
The  availaUe  data  arid  the  rpptoximations  which  worked  so  nicely  die  aUudies  [2-4]  are  insufficiendy 

accurate  for  CO.  In  particular,  die  atomic  electric  transition  dqxde  moments  (whether  obtained  fiom  theory 
or  atperinient)  are  highly  uncertain.  We  have  noted  that  even  wdl-conelated  computed  transition 
moments  P8]  may  disagree  significantly  with  experimental  results.  For  this  mote^e,  which  is  probably 
a  "worst  case"  for  the  BAF  qiproach,  odwr  wrys  to  rpproximate  die  dispersion  energy  may  be  more 
reliaUe.  Clearly  dien,  we  have  not  obtained  a  truly  satisfactory  descr^tion  of  the  di^rsion  interaction 
energy.  However,  in  a  systematic  study  of  rare  gas  interactions  Krauss  and  coworkers  [1]  have  found  that 
the  neglect  of  overlq)  effects  leads  to  an  overestimation  of  die  dispersion  mergy.  Thus,  our  neglecting 
overlap  effects  in  the  induced  dqrole-induced  dqxiie  interaction,  v^hich  we  do  tr^  conqiensates,  at  least  in 
part,  for  our  neglecting  die  higha  order  terms. 

With  the  possible  exception  of  the  1  ^2'*’  state,  the  (A,S)  coiqiling  case  is  not  accurate  enough  even 
at  R  =  7  a.u.  due  to  the  strong  spin-orbit  interaction  in  the  oxygen  atom.  So,  spin-orbit  interactions 
must  be  induded  and  die  (JJ)  coupling  case  must  be  considered  for  an  accurate  description  of  the  long- 
range.  In  that  coupling  scheme,  long-range  bound  states  are  predicted  for  some  "1"  ar^  "2"  symmetry 
mdecular  states.  Unfortunately,  no  experimental  data  to  connate  these  with  are  known  for  these  long- 
range  molecular  states. 
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Table  1 :  Mean-square  radii  <r2>/ao2  of  incomplete  p  shells  for  carbon  and  oxygen  3p 
term  and  R/ao^  resulting  value  (ao  is  the  Bohr  radius) 


Ref 

Carbon 

Oxygen 

R 

This  woik 

3.74717 

1.97801 

8.89433 

H.F12 

3.746801 

1.974975 

8.8798 

H.F12 

3.889709 

2.001392 

9.3418 

H.F.S13 

3.892540 

1.93133 

9.0213 

D.Fl^ 

3.88985 

2.00067 

9.3387 

H.Fll 

3.890 

2.001 

9.341 

Re£-7 

3.7483 

1.9741 

8.8794 
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Table  2a:  A  comparison  of  line  strengths  for  some  electric  dipole  transitions  in  carbon 


Transition 

Present 

F.F.I39] 

Other  theories 

Expt. 

^expt(  A) 

2s2p3  3D0 

5.21 

4.7lf.  3.75*.  3.56® 

4.2h.  3.84a-b.  3.52c-<l 

1561.1 

28^3po 

— 

12.47 

3.82f*  3.62* 

1.5‘>.  2.484 

1329.3 

2s^3so 

4.446 

7.68.  7.62* 

945.5 

(2po)3s  3po 

6.74 

— 

5.30f 

8.3h.  6.52C 

1657.2 

(2po)4s  3po 

0.84 

— 

— 

0.76*»  ,1.864 

1280.4 

(2po)3d3po 

0.65 

l.lh.  1.42» 

1261.3 

(2po)3d3DO 

2.3 

— 

— 

2.4h.  3.406i 

1277.5 

a-D.  T.  Pegg  et  al.  (1970)  [15] 
b-M.  C.  Poulizac  et  al.  (1971)  [16] 
c-G.  M.  Lawrence  et  al.  (1966)  [17] 
d-N.  H.  Brooks  et  al.  (1977)  [18] 
e-C.  A.  Nicolaides  et  al.  (1971-73)  [21,34] 
f-A.  W.  Weiss  (1967)  [22] 
g-A.  B.  Bolotin  et  al.  (1956)  [20] 
h-O.Boldt(1963)  [19] 
i-J.  Bronuukler  et  al.  (1978)  [40] 
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Table  2b:  A  comparison  of  line  strengths  for  some  electric  dipole  transitions  in  oxygen 


Transition 

Present 

F,F(2S,32) 

Other  Theories 

Expt. 

^exp( A ) 

2s2p5  3po 

5.30 

4.64S.  1.63^  4.46't. 

6.69^.® 

792.0 

2MK  3.45*> 

(^SO)3s  3S0 

1.26 

1.75 

0.91«.  ZIO**,  1.2CJ 

1.85“.  1.81®.* 

1303.5 

4s 

0.22 

— 

0.16®,  0.29^ 

0.42®,  0.28* 

1040.1 

5s 

0.081 

0.057®.  0.0974 

0.23® 

977.2 

6s 

0.040 

0.0027®.  0.0464 

— 

951.6 

7s 

0.022 

— 

0.015® 

— 

938.5 

8s 

0.014 

— 

9.28E-3® 

930.9 

9s 

9.35E-3 

— — 

6.11E-3® 

— 

926.1 

10s 

6.57E-3 

— 

4.24E-3® 

922.9 

(2d«>)3s  3i>> 

2.4 

1.49 

1.654,  1.6i.® 

1.7“,  1.43®,  1.49*  989.5 

(2po)3s  3po 

1.62 

1.01 

2.064,  0.96i.c.  1.92*>,  2.13“, 

3.27“,  4.0lf 

878.5 

(^»)3d  3d® 

0.33 

— 

0.33®,  0.634,  0.3 li 

0.88® 

1026.6 

4d 

0.16 

— 

0.434, 0.17® 

0.34® 

972.5 

5d 

0.09 

— 

0.092®,  0.184 

0.21® 

949.4 

6d 

0.053 

— — 

0.054®,  0.114 

0.15® 

937.3 

7d 

0.033 

0.035® 

0.12® 

930.2 

8d 

0.022 

— — 

0.023® 

- - 

925.6 

9d 

0.016 

0.017® 

922.5 

(2D®)3d3po 

0.22 

— 

0.0374,  0.18j,  0.33“  0.86“ 

6.4lf 

811.4 

a-N.  H.  Brooks  et  al.  (1977)  [18] 

b-W.  H.  Smith  et  al.  (1971)[23] 

c-B.  H.  Armstrong  (1966)  [26] 

d-A.  K.  Pradhan  and  H.  E.  Saraph  (1977)  [27] 

e-S.  Chung  et  al  .(1986)  [28] 

f-E.  J.  Knystautas  et  al.  (1973)  [29] 

g-U.  I.  Saftonova  (1975)  [33] 

h-W.  L.  Luken  et  al.  (1976)  [36] 

i-C.  Zeippen  et  al.  (1977)  [37] 

j-P.  S.  Kelly  (1964)  [24] 

k-C.  A.  Nicolaides  (1973)  [34] 

1-C.  A.  Nicolaides  et  al.  (1976)  [41] 
m-O.  M.  Lawrence  (1970)  [30] 
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Table  3:  Long-range  coefficients  C5  and  C6  for  the  interactions  C(3p)  +  0(3p) 


Molecular  State 

Footnote 

C5  (a.u.) 

C6  (a.u. 

1  1,342+ 

a 

-11.20 

-5.81 

b 

-11.20 

-7.34 

c 

-11.20 

-8.76 

1  U4a 

a 

-1.87 

-5.16 

b 

-1.87 

-7.47 

c 

-1.87 

-8.51 

2  1.342+ 

a 

0. 

-4.85 

b 

0. 

-6.42 

c 

0. 

-7.28 

1  14.32- 

a 

0. 

-5.61 

b 

0. 

-7.92 

c 

0. 

-9.29 

1  144n 

a 

0. 

-4.81 

b 

0. 

-6.08 

c 

0. 

-6.97 

2l44n 

a 

+7.47 

-5.70 

b 

+7.47 

-6.95 

c 

+7.47 

-8.49 

^Extended  hydiogenic  approximation  for  all  dipole  transition  moments  except: 

a.  G.  A.  Nicolaides  values  [34]  for  the  2s^  2p°  -  2s  2p°^^  transitions. 

b.  G.  Froese-Fischer's  values  [32,29]  for  the  2s22p°  -  2s  2p“+*  transitions. 

c.  Correlated  calculations,  experimental  or  estimated  values  for  the  2s^2p^  -  2s^2p°*l  3s,  3d  transitions 
and  values  given  by  G.  Froese-Fischer  for  the  28^  2p“  -  2s  2p°+^  transitions. 
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Table  4:  Position  and  Interaction  Energy  of  the  long-range  extrema  of  CO 
dissociating  in  C(3pj)  +  0(3pj) 


Molecular 

Interaction 

Symmetry 

State  no 

Extrema  Type 

Position  (a.u.) 

Energy  (cm*^) 

0+ 

7 

max 

~7. 

+3.4 

8 

mm 

~8.8 

-2.7 

max 

-«.8 

+1.2 

0- 

2 

max 

-9.7 

+1.9 

9 

min 

-8.2 

-2.3 

1 

7 

max 

—6.6 

+24.9 

9 

min 

-7.6 

-11.8 

max 

-7.0 

-11.1 

10 

max 

-6.6 

+15.4 

11 

min 

-7.0 

-6.6 

max 

-6.4 

-3.2 

12 

min 

—8.6 

-1.7 

2 

5 

max 

-8.6 

+3.1 

6 

max 

-6.6 

+29.8 

8 

min 

-7.7 

-15.5 

max 

-6.4 

-4.0 

10 

min 

-8.1 

-2.7 

3 

3 

max 

-6.5 

+33.6 
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Appendix  I:  Molecular  basis  functions  for  C(%) -t- 0(^P) 
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31*  =~  [1110(^.^11 14(bj,^114(b.^lllOO>^ 
51*  =llllO>Jlll(fc> 
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^[111 10>J1  l-l±l>^j-ll  100>^.ll  10±l>j,+ll  1-10>^U  1 1±1>^ 


II 1 10x11 10±lx-ll  100>J1 1 1±1> 


6 1+  2x(»3^I+) 


Appendix  II:  Molecular  basis  fiinctions  for  C(3pj)  +  0(3pj)  interaction 


1 .  ^Po(C)  +  ^P2(0)  interaction:  Edhs  =  0.0  cm-' 
0+  =  I005cI20^ 

I=±l0(bcl2±l50 

2*I00>cI2±25o 

2.  ^Pi(C)  +  ^PilO)  interaction:  Ediss  =  16.4  cm*l 

0+  =  I  I  l-l^c  1 2150- 1 11>C  '  2-150  ) 

0-  =  ;^  {  I  1-15C  I  21^0  +  I  1 1^  I  2-150  ) 

O’  =  I  lOio  *  2050 
1  =  1 1  ±  15C  I  2050 

1  =11050 1 2±  150 
1=1  1+15cI2±25o 

2  =  ±l  l±l5cl2±l5o 

2  =  ±I105cI2±25o 

3  =  1 1  ±  150  I  2  ±  250 

3.  ^P2(C)  +  ^PilO)  interaction:  Edh*  =  43.4  cni*l 
0^'  =  ;^|f2 -250  I  2  250  +  1 225012-250) 

0+  =  ;^!  12- 15012  150  +  12150  I  2- 150) 

0-  =  ^{I2-25oI225o-I225oI2-25o) 

0-  =  ;^{l2-l5ol2  1:o-'2l5cl2-l5o) 

0+  =  I205oI205o 
I=±l2±l5ol2  05o 
1  =±120501  2±  l5o 
1  =  ±  1 2+150  I  2  ±  250 


3  states 


9  states 


14  states 
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2  =  I2±25cI  20>o 
2  =  I2±1^  l2±l>o 


2  =  I20?cI2±23o 

3  =  ±l2±25cl2±l5o 

3  =  ±I2±1:cI2±25o 

4  =  I2±25cI2±2>o 

4.  ^Po(C)  +  ^Pi(O)  interaction:  Ediu  =  158.3  cni"^ 
0-  =  l005cl  1  050 

l  =  l00^ll±l>O 

5 .  ^Pi(C)  +  ^Pi(O)  s  Edijg  interaction:  174.7  cni*l 
0+  =  ^{  11-1^1  1  130  +  111^11-150) 

0+  =  1 103c  I  1050 

0-  =  ~(ll-l3cl  1  l:o-lll5Cll-l=^l 
l=±ll±l3ci  1  050 

1  =  ±l  1  03c  I  1  ±  l50 

2  =  1 1  ±  13C  I  1  ±  1>0 


2  states 


6  states 
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6.  ^P2(C)  +  ^Pi(O)  interaction:  Ediss  =  201.7  9  states 


0+  =  ;~{l21:cll-l5o-l2-l>cll  1=0) 

0-  =  ;^|I21>cI1-1>O  +  I2-1>cI1 

0’*l20>c  I  1  03q 

l=l20:cl  1±1>0 
l=l2±l>cl  1  050 

1  =  I2±25cI  1+1>o 

2  =  ±l2±l5cl  1±  l5o 
2  =  ±I2±25cI  1  050 
3*I2±25cI  1  ±1^ 

7.  3Po(C)  +  3Po(0)  interaction  =  EdUs  =  227  cm-l  1 

0+  =  l005c  10  050 

8.  3pj(C)  +  3Po(0)  interaction  a  Ediw  =  243.4  cm-l  2 

0-sll05ol0  0:O 

1  =  1 1  ±  l5c  I  0  050 

9.  ^PilC)  +  3p,(0)  interactions  =  Edtss  r  270.4  cm’*  3 

0+  =  1 2050  *0  050 
la±l2±  15010  050 

2  =  1 2±2:o'  0  050 


Total:  43  molecular  states 


state 


states 


states 
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Appendix  III:  2p  Orbital  SCF  Wavefunctions  for  Carbon  and  Oxygen 
Atoms 

'F{2p)  =  ^  Ni  Q  r  Ni-l  e-  Cj  Y™*  (e,d>) 
i  ' 


Atom 

n 

1 

Ci(2pa) 

C|  (2pa) 

C 

2 

1 

6.51 

-.010247216 

-.0102472161 

2 

1 

2.6005 

-.2319540338 

-.2319540339 

2 

1 

1.4436 

-.5469976538 

-.5469976544 

2 

1 

.9807 

-.2823784197 

-.2823784190 

O 

2 

1 

7.5648 

-.0765520772 

-.0165527165 

2 

1 

3.4499 

-.3277071772 

-.3225525816 

2 

1 

1.8173 

-.6188261938 

-.5656618839 

2 

1 

1.1439 

-.1184586922 

-.1867904703 

2s  Orbital  HF  Wavefunctions  for  Carbon  and  Oxygen  Atoms 

Atom 

n 

I 

Cl 

Cl 

C 

1 

0 

9.4826 

.0107111691 

1 

0 

5.4360 

.2081422862 

2 

0 

4.2010 

.1476498490 

2 

0 

2.6844 

-.3354681690 

2 

0 

1.5243 

-.7504500909 

2 

0 

1.0575 

-.0810066997 

0 

1 

0 

13.7574 

-.0047624085 

1 

0 

7.6141 

-.2215772275 

2 

0 

5.8660 

-.1918271486 

2 

0 

4.3120 

.2536208043 

2 

0 

2.4802 

.6082459022 

2 

0 

1.6982 

.3482806491 
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Appendix  IV :  Explicit  E5  and  E5  Matrix  Elements  for  and  ^3.3^  molecular  states  in 

tennsof9l^and5lap  witha,  PE  {0,1^}  9l^and5l<xp  *91^  are  given  by  the 
relations  (1)  and  (2)  respectively. 


1.  1.3,52+  states 

rE5(l.l)  E5(1.2)1  .5656854  ]  0,2 

E5(2.2)  _  _  T(5  _ 


E5(2,2) 
E6(1.1)  E6(1,2) 
E6(2.2) 


I  .31426968  ^ 


-.0628539  _ 

8  (91o2  +  9^20) 


1^  .0785674 
+  ^  5  9ln  + 


2.  1,3,52-  states: 


4714045  _  ™  .040855 

1  (9li2  +  9l2l)+  ^  41  9 

“3?y  J  L  ^  . 


E3(1,1)*C5(1)*=0 

E6(1,1)  =  C3(l)  =  -  ^(9loi  +  9lio)  +  ^9lii  -  ^(91 12  +  9I21)  +  7^9122 


3.  i,3,5n  states: 

-Esd.l)  E3(1.2)- 

E5(2,2)  . 
"E6(1,1)  E6(1,2)“ 
E6(2,2)  . 


4  4  1 

■nr  1?!  L,2«.2 


4 


5  11  r  I  2  " 

2  9ioi  +  4  "02  + 

?  J  L  5?T  . 


■  25 

55jy 

0 

25 

ml 

9111  + 

95  2 

”  5WT 

35 


9112+ 


4  2 

1 

31T 


9^20  + 


2  1  1 
9  7 

5  5lio 

m  J 

35  2  1 

-m  ^  gf  ^ 
95  ^21  + 

'  mJ 


4 ,  1,3, 5 A  states: 


E3(l,l)-C5(l)-^9tj9tb 


E«(l,l)  *  C«(l)  »  ^(3*02  +  3*20)  +  ff  3i|i  -  w(91i2  +  3I21)  +  ot3*22 


